The Role of Reconstructed Surfaces in the Intrinsic Dissipative Dynamics of Sihcon 

Nanoresonators 
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Dissipation in the flexural dynamics of doubly clamped nanomechanical bar resonators is inves- 
tigated using molecular dynamics simulation. The dependence of the quality factor (Q) on temper- 
ature and the size of the resonator is calculated from direct simulation of the oscillation of a series 
of Si (001) bars with bare {100} dimerized surfaces. The bar widths range from 3.3 to 8.7 nm, all 
with a fixed length of 22 nm. The fundamental mode frequencies range from 40 to 90 GHz and Q 
from 10^ near 1000 K to lO'' near 50 K. The quality factor is shown to be limited by defects in the 
reconstructed surface. 
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There is currently much interest in ultrahigh frequency 
nanoscale mechanical resonators for both fundamental 
science and engineering applications. From a fundamen- 
tal perspective, nanoelectromechanical systems (NEMS) 
are leading candidates for observing macroscopic quan- 
tum behavior P, Q , while applications include force sens- 
ing and frequency generation |3| • Achieving the quan- 
tum limit requires resonators with fundamental mode 
frequencies uj in the radio-to-microwave frequency range 
that are also well separated from other low-lying modes 
on the scale of the resonance frequency width. The fre- 
quency width is determined by the mechanical energy 
damping rate 7 and hence large quality factors Q = uj/j 
are required. Force sensing and frequency generation ap- 
plications also require large quality factors. For exam- 
ple, the relative uncertainty in a generated frequency is 
1/Q [5]. It has also recently been demonstrated that 
nanoresonators can be used as mass detectors with sensi- 
tivities approaching the single Dalton (i.e., carbon atom 
mass) level [1, 0, U0] ; when a mass attaches to the res- 
onator, the fundamental resonance frequency shifts and 
the ability to resolve this shift depends on the quality 
factor. 

The energy dissipation of nanomechanical resonators 
is described in general by three mechanisms [l^, [lH: (1) 
energy losses to the extrinsic environment through the at- 
tachment cross section of a resonator to its support sub- 
strate ('clamping loss'), or through air friction; (2) intrin- 
sic losses arising from anharmonicity (or phonon scatter- 
ing), thermoelastic damping, electron-phonon coupling 
and the strain-driven motion of defects in the interior 
(bulk) of a resonator; and (3) intrinsic losses arising from 
analogous processes within the surface layers. For exam- 
ple, lattice defects can move in response to an applied 
strain and then relax via elastic radiation. The energy 
dissipated depends on both the microscopic structure of 
the defect dynamics and on the long-wavelength acoustic 
properties of the mechanical device and its surroundings. 
The quality factors of the fundamental flexural and other 
low-lying vibrational modes have been observed to de- 



crease approximately linearly with decreasing resonator 
size [3]. This size dependence is usually interpreted to 
be a signature of surface-related dissipation mechanisms 
being dominant in (sub-)micron scale resonators. 

In this letter, we investigate the fundamental, intrin- 
sic energy dissipation for doubly clamped silicon bar 
nanoresonators arising from the anharmonic interatomic 
potential and surface reconstruction; the resonators are 
in vacuum and support atoms kept fixed so that extrinsic 
loss mechanisms are suppressed. Direct numerical simu- 
lations of the dynamics of individual atoms are employed 
to examine how energy is transferred from the fundamen- 
tal mode of interest to the other modes. We have exten- 
sively adapted the scalable, parallel molecular dynamics 
(MD) code MDCASK [13] to simulate and analyze the 
nanoresonator oscillation, making use of the Stillin ger- 



Weber (SW) potential [13| and the Tersoff potential [14 1 
to describe the interatomic interaction. These potentials 
involve nearest-neighbor two-body and three-body inter- 
actions parameterized empirically based on bulk elastic 
and thermodynamic properties of the Si diamond-cubic 
crystal and, for SW, its melt point. Both potentials are 
widely used and validated, achieving a reasonable de- 
scription of many properties including surface structures 
and energies [3] ■ The use of both potentials in our sim- 
ulations enables their systematic comparison in the mod- 
eling of mechanical dissipation in Si nanoresonators, for 
which their relative merits are not known a priori. Sev- 
eral authors ITtI . [ist have used MD simulations to 



study energy dissipation in carbon nanotube oscillators. 
Our MD study of Si resonators is conducted in a man- 
ner similar to that used previously for the oscillations of 
quartz resonators l^, 20j . 

The simulations focus on single-crystal diamond-cubic 
Si (001) bars with {100} faces. To accommodate the sym- 
metric dimer 2x1 surface reconstruction, square cross 
sections of m x m unit cells are selected with even m: 
m = 6, 8, . . . , 16, corresponding to widths of 3.26, 4.34, 
5.43, 6.52, 7.6 and 8.69 nm. The bar length is fixed to 
be 40 unit cells (~ 22 nm). Including the dimer atoms 
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FIG. 1: (Color) Time evolution of the center of mass in 
the transverse x (blue), y (red), and longitudinal z (green) 
directions (in Angstroms) for the 8.69-nm Tersoff bar near 
1000 K. The black solid curve is the root-mean-square ampli- 
tude whose log is used to extract the damping rate. 
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FIG. 2: (Color) The left half of the 5.43-nm Tersoff bar near 
1000 K after 3000 ps. Atoms are colored according to their 
potential energy. 



on the surface, the number of atoms in these bars ranges 
from 12,360 to 84,360 atoms. The bars are initially ther- 
malized and mechanically equilibrated at 50 K, 100 K, 
200 K, 300 K and 1000 K. The thermostat is then turned 
off for the duration of the simulation: 3,000 ps (3 million 
time steps). In order to clamp both ends of the bars, 
the positions of the atoms initially within one unit cell at 
both ends of the bar are kept fixed. The lowest flexural 
mode is excited by adding an initial transverse velocity 
in the x-direction to each atom on the resonator, varying 
with z according to the normal mode derived from linear 
continuum elasticity [21]. The initial amplitude is chosen 
to be 10% of the bar width. 

When the bars are "plucked" in the x direction, the 
x-mode energy is transferred partially over time into the 
y-direction mode, as evidenced by the center of mass tra- 
jectory shown in Fig.[TJ The energy transfer between the 
two transverse modes is attributed to the fact that the 
silicon dimer atoms on the bar surface {100} align them- 
selves in the diagonal direction as shown in Fig. [21 The 
dimerized surface of the bar breaks the cubic symmetry 
of the bulk and, we infer, induces the mixing of the two 
transverse modes. A tiny amount of the transverse mode 
energy is also transfered into the longitudinal mode in 
the z direction, but the transfer is very limited due to 
the clamping at both ends. 

The quality factor Q is calculated as the ratio of the 
mode frequency to the energy damping rate, uj/j, where 
the frequencies of the dominant amplitude modes are 
identified from the peaks of the Fourier transform of the 
center of mass position. The frequencies for the two dom- 
inant transverse modes are very close to each other, giv- 
ing rise to the beat pattern observed in Fig.[TJ The damp- 
ing rate is extracted as twice the slope in the time trajec- 
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FIG. 3: (Color) Log of the amplitude A{t) relative to the 
initial amplitude ^(0) versus time t at various temperatures 
for the 6.52-nm Tersoff bar. The dotted curves are from the 
simulations and the straight solid lines in blue are the fits 
whose slopes correspond to one half the energy damping rates. 



tory of the log of the envelope of the amplitude which is 
found to decay exponentially at a constant rate to a good 
approximation, as shown in Fig.[3l Due to the mode mix- 
ing, the amplitude is taken to be the root-mean-square of 
the sum of the mean-adjusted center-of-mass amplitudes 
in all three directions. 

The temperature dependence of the damping rates for 
the various bars is plotted in Fig. [H The fundamental 
mode frequencies range from 40 to 90 GHz, increasing 
with bar thickness and decreasing slightly with temper- 
ature. As the temperature increases, the quality factor 
decreases from lO'^ near 50 K down to 10^ near 1000 K. 
The room temperature values are of the order of, but 
somewhat larger than, the highest quality experimental 
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FIG. 4: (Color) Temperature dependence of the damping 
rates, plotted in terms of ln(Q/a;) versus l/T for Tersoff bars 
in blue and Stillinger- Weber bars in red. uj is the mode 
frequency in radians/ps. The curves are the best fits for 
ln(7i +72)"^. The temperature of each data point represents 
an average over the whole simulation period. The error bars 
are estimated from the uncertainty in the slope coefficient for 
the least-squares fitting of a straight line. 



microrcsonators scaled to the simulated sizes Q. Quan- 
titatively, the damping can be attributed to an ather- 
mal process which dominates at low temperatures and a 
thermal relaxation process which is manifested at room 
temperatures. Denote the corresponding damping rates 
for the two processes to be 71 and 72 respectively: 



QiT) 



coiT) 



71 



72 (T)' 



(1) 



where u> is the mode frequency of the resonator and the 
athermal process is assumed to have the constant damp- 
ing rate 71 for both Tersoff and Stillinger- Weber bars. 
72 describes the thermal relaxation process in a standard 
linear solid and behaves according to Arrhenius' law [2^ . 



72(T) «7oexp(~i?/fcBr). 



(2) 



where H is the activation energy and 70 is the constant 
damping rate in the high temperature limit. 

The curves in Fig. 2] give the best fits of this functional 
form to the two data sets. The fitted athermal damping 
rate is 71 w 2.75 x 10^^ ps^^ for both interatomic po- 
tentials. The thermal damping rate is estimated to have 
70 w 1.83 X 10-3 ps-i and H « 50.9 meV for the Tersoff 
case while 70 « 9.22 x 10"^ ps"^ and H sa 57.9 meV 
for the Stillinger- Weber case. When extracting these fit- 
ting parameters we have excluded the Stillinger- Weber 
bars near lOOOK where an additional relaxation process 
appears to take place. We will discuss the high temper- 
ature regime in more detail below. 

Over the course of the simulation various degrees of 
surface roughening are observed as shown in Fig. [21 In 
most bars the roughening on the surfaces and the edges is 



FIG. 5: (Color) Potential energy of randomly selected sample 
atoms in the 5.43-nm Tersoff bar at ~1000 K over the first 
100 ps. The three energy bands, denoted as the edge (blue), 
the dimer (red) and the bulk (cyan) , each comprise 20 sample 
atoms. 



more pronounced at higher temperatures and there is less 
roughening in the Tersoff bars, which correspondingly 
have less damping and higher Q's than the Stillinger- 
Weber bars. At higher temperatures, a greater portion 
of the edge region was roughened as compared with the 
interior surface and bulk regions. This observation is 
consistent with theoretical and experimental studies on 
edge roughening of equilibrium crystal shapes, which can 
be explained by the thermal relaxation of the steps on 
the surface [23!, [13, ■ The activation energies extrapo- 
lated from Fig. m are smaller than the activation energies 
for various diffusion processes of 2 x 1 Si dimers on the 
{100} surface, ranging between 0.7 eV to 1.36 eV [2^ . 
In fact they are comparable to step creation energies; 
e.g., 10 meV/a for the single-layer step, denoted as Sa, 
or 50 meV/a for the double-layer step, denoted as Db, 
where type A{B) refers to the steps with dimer rows on 
the upper terrace oriented parallel (perpendicular) to the 
step edge [27| and a « 3.85 A is the 1x1 surface lattice 
constant. Hence the dissipation is very likely dominated 
by the thermal relaxation of the step edges and not by 
the diffusion of dimers. 

The coloring of the atoms according to their potential 
energies in Fig. [2] suggests that the atoms cluster into 
three energy "bands" , denoted here as the edge atoms 
whose energy are greater than -3 eV, the bulk atoms 
with energy less than -4 eV and the surface dimer atoms 
with energy between -3 eV and -4 eV. Figure [5] is a sam- 
ple plot of randomly selected atoms in each band. Most 
atoms remain within the same energy band throughout 
the simulation. As temperature rises, more atoms hop 
between bands. All the atoms that undergo the largest 
energy hopping are located within one unit cell or so 
from the geometric edges. The work done by the funda- 
mental mode on the atoms during these hops is released 
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irreversibly and leads to dissipation. We find that the 
activation energies associated with the hopping rates are 
of the same order of magnitude as those extracted from 
Fig.H 

Another observation is that the fitted damping rates 
in Fig. m do not seem to depend on the bar size. This 
invariance is another indication of the quality factor be- 
ing limited by defects near the geometric edge regions. 
When the edge defects dominate the dissipation, differ- 
ent cross section but equal length bars are expected to 
have similar dampin g ra tes at a give n temperature. Cur- 
rent experiments (sl. l28l. 29l . 30l . 31|, on the other hand, 
involve larger silicon resonators than those in our simu- 
lations. These resonators, usually made from a top-down 
approach, are expected to have vicinal surfaces with step 
edges in addition to the geometrical edges. The typical 
width between reg ularly spaced step arrays ranges be- 



tween 2-10 nm [3^|33|. When the resonator dimensions 



are large as compared with the step spacings, the damp- 
ing rate is expected to depend on the surface area. 

At room temperature the vicinal silicon surface shows 
a double-step reconstruction above a critical miscut angle 
of 3° while at smaller miscut angles the single-layer steps 
dominate the surface 13, 35 1. For surface with a large 
miscut angle, a phase transition from Db to S*^ steps has 
been observed near 800 K [3^ . Hence the smaller forma- 
tion energy of Sa steps may be related to the smaller 
slope in Fig.|4]for the Stillinger- Weber bars near lOOOK. 

In conclusion, we have conducted the first systematic 
study of dissipation and the resulting quality factors of 
nanoscale bar resonators using molecular dynamics simu- 
lations, finding that the quality factors of the Si nanores- 
onators are dominated by thermal relaxation of the step 
edges. Furthermore, transfer of energy between modes 
with polarization transverse to the {100} surfaces is ob- 
served, due to symmetry breaking of the bulk by the 
dimer surface reconstruction. The present results further 
reinforce the relevance of the surface for the mechanical 
properties of nanoscale resonators. 
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